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Abstract. 

We apply the perturbative chiral quark model (PCQM) at one loop to analyze 
the axial form factor of the nucleon. This chiral quark model is based on an 
effective Lagrangian, where baryons are described by relativistic valence quarks and 
a perturbative cloud of Goldstone bosons as dictated by chiral symmetry. We apply 
the formalism to obtain analytical expressions for the axial form factor of the nucleon, 
which is given in terms of fundamental parameters of low-energy pion-nucleon physics 
(weak pion decay constant, strong pion-nucleon form factor) and of only one model 
parameter (radius of the nucleonic three-quark core). 
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1. Introduction 

The nucleon axial form factor is of fundamental significance to weak interaction 
properties and to the pion-nucleon interaction. Hence it provides one important 
test for theories that attempt to describe the structure of the nucleon. The current 
status of the experimental and theoretical understanding of the axial form factor of 
the nucleon is reviewed in Based on the original development of chiral quark 

models with a perturbative treatment of the pion cloud [H]-[E], we recently extended 
the relativistic quark model suggested in |S] to describe the low-energy properties of 
the nucleon In references |H]-[IE| we developed the so-called perturbative chiral 

quark model (PCQM) in application to baryon properties such as: sigma-term physics, 
electromagnetic form factors of the baryon octet, nN scattering and electromagnetic 
corrections, strange nucleon form factors, electromagnetic nucleon-delta transition, etc. 
In the present work we follow up on the earlier investigations and employ the same 
model in order to study the axial form factor of the nucleon. 

The PCQM is based on an effective chiral Lagrangian describing quarks as 
relativistic fermions moving in a self-consistent field (static potential) V e g(r) = S(r) + 
7°V(r) which is described by a sum of a scalar potential S(r) providing confinement 
and the time component of a vector potential 7°l / (r). Obviously, other possible Lorenz 
structures (e.g., pseudoscalar or axial) are excluded by symmetry principles. It is known 
from lattice simulations that a scalar potential should be a linearly rising one and a 
vector potential thought to be responsible for short-range fluctuations of the gluon field 
configurations [16 . In our study we approximate V e s(r) by a relativistic harmonic 
oscillator potential with a quadratic radial dependence jU] 



The model potential defines unperturbed wave functions for the quarks, which 
are subsequently used to calculate baryon properties. This potential has no di- 
rect connection to the underlying physical picture and is thought to serve as an 
approximation of a realistic potential. The vector part of the potential is also a pure long- 
ranged potential and is not responsible for the short-range fluctuations of gluon fields. In 
general, we need a vector potential to distinguish between quark and antiquark solutions 
of the Dirac equation with an effective potential. Note, that this type of the potential 
was extensively used in chiral potential models |H],[IZ|-[221- A positive feature of this 
potential is that most of the calculations can be done analytically. As was shown in 
Refs. |H],[IZ|-[221 and later on also checked in the PCQM [Oj-jlSl, this effective potential 
gives a reasonable description of low-energy baryon properties and can be treated as a 
phenomenological approximation of the long-ranged potential dictated by QCD. 

Baryons in the PCQM are described as bound states of valence quarks supplemented 
by a cloud of Goldstone bosons (n, K, rj) as required by chiral symmetry. The chiral 
symmetry constraints will in general introduce a nonlinear meson-quark interaction, but 
when considering meson as small fluctuations we restrict the interaction Lagrangian up 
to the quadratic term in the meson fields. With the derived interaction Lagrangian 





(1) 
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we do our perturbation theory in the expansion parameter 1/F (where F is the pion 
leptonic decay constant in the chiral limit). We also treat the mass term of the current 
quarks (m, m s ) as a perturbation. Dressing the baryonic three-quark core by a cloud 
of Goldstone mesons corresponds effectively to the inclusion of sea-quark contributions. 
All calculations are performed at one loop or at order of accuracy 0(1/F 2 , m, m s ). 

To be consistent, we use the unified Dirac equation with a fixed static potential both 
for the ground and for the excited quark states. In the Appendix we give details of the 
solutions to the Dirac equation for any excited state. Inclusion of excited states should be 
handled consistently. First of all, one should guarantee conservation of local symmetries 
(like gauge invariance). Second, excited states should be restricted to energies smaller 
than the typical application scale A w 1 GeV of low-energy approaches. An alternative 
possibility to suppress the inclusion of higher-order excited states is to introduce a 
meson-quark vertex form factor [HI EH- Solving the Dirac equation with a relativistic 
harmonic oscillator potential ((TJ one can show, that the energy shift between the first 
low-lying lp excited states and the Is ground state is about 200 MeV. The excited 
states (Id and 2s) lie about ~ 370 MeV above the ground state. The 2p and the 1/ 
states are 530 MeV heavier when compared to the ground state. As soon as the typical 
energy of the ground-state quark is about 540 MeV^one can restrict to the low-lying 
lp, Id and 2s excited states with energies smaller than the typical scale A = 1 GeV. 
The requirement of convergence of physical observables when including excited states 
is physically not meaningful since it takes states with very large energies where the 
phenomenological low-energy approaches break down. One of us showed in Ref. jBJ that 
the inclusion of excited states to the nucleon and A masses can be convergent when using 
a linearly rising confinement potential, e.g., the use of potential with a quadratic radial 
dependence leads to a nonsatisfactory convergence. However, this statement is sensitive 
to the quantity one is testing. On the other hand, our approach has some different 
features in comparison to previous ones (see, e.g., Ref. [E!,[II]-!22])- ^ n particular we 
perform a consistent renormalization procedure when we include meson cloud effects. It 
gives additional contributions to physical quantities which were not taken into account 
before. In Ref. [H] we demonstrated that excited quark states (lp, Id and 2s) can 
increase the contribution of loop diagrams but in comparison to the leading order (three- 
quark core) diagram this effect was of the order of 10%. This is why we were interested 
to study these effects for the example of the axial form factor. Again, we truncate the 
set of excited states to the lp, Id and 2s state with energies which satisfy the condition 
£ < A = 1 GeV. We do not pretend that we have a more accurate estimate of the 
whole tower of excited states. The scheme we use is thought to take the excited states 
into account in an average fashion. We showed that the zeroth-order value of the axial 
nucleon charge is not changed much in the presence of meson cloud effects in consistency 
with chiral perturbation theory. As an extention of previous work, which is dominantly 
aimed to describe the low-energy static properties of baryons, we consider the model 

f This value can be deduced from a calculation of octet and decouplet baryon spectrum. Similar 
estimates can be found in other chiral quark calculations. 
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prediction for the axial charge and for completeness also for the axial form factor. No 
further parameters are adjusted in the present work. 

In the present paper we proceed as follows. First, we describe the basic features 
of our approach: the underlying effective Lagrangian, the unperturbed, that is valence 
quark, result for the nucleon description together with the choice of parameters and a 
brief overview of perturbation theory when including the meson fields. The full details of 
the renormalization technique can be found in reference jU] . In section III we concentrate 
on the detailed analysis of the nucleon axial form factor in our approach. We derive 
analytical expressions in terms of fundamental parameters of low-energy pion-nucleon 
physics (weak pion decay constant, strong pion-nucleon form factor) and of only one 
model parameter (radius of the three-quark core of the nucleon). Numerical results in 
comparison with data are presented to test the phenomenological implications of the 
model. Finally, section IV contains a summary of our major conclusions. 



2. The perturbative chiral quark model 

Following considerations lay out the basic notions of the perturbative chiral quark 
model (PCQM), a relativistic quark model suggested in Ej and extended in 

for the study of low-energy properties of baryons. In this model quarks move in 
an effective static field, represented by a scalar S(r) and vector V(r) component with 
^eff(r) = S(r) + 7°y(r) and r = \x\, providing phenomenological confinement. The 
interaction of quarks with Goldstone bosons is introduced on the basis of the nonlinear 
cr-model [2Bj- The PCQM is then defined by the effective, chirally invariant Lagrangian 

Anv maun 

C im {x)=$(x){i?-<fV{r)-S(r) 

+^Tr[d^Ud^], (2) 

with an additional mass term for quarks and mesons 

C x sb(x) = -${x)Mil>{x) -jTt[^ 2 M}, (3) 

which explicitly breaks chiral symmetry. Here ip is the quark field; U = e 1$//F is the 
chiral field; $ is the matrix of pseudoscalar mesons (in the following we restrict to the 
SU{2) flavor case, that is $ —* n = iz ■ r); F = 88 MeV is the pion decay constant in 
the chiral limit; M. = diag{m,m} is the mass matrix of current quarks (we restrict to 
the isospin symmetry limit with m u = ma = rh = 7 MeV) and B = — {0\uu\0)/F 2 = 1.4 
GeV is the quark condensate parameter. We rely on the standard picture of chiral 
symmetry breaking and for the mass of pions we use the leading term in their chiral 
expansion (i.e. linear in the current quark mass): M% = 2rhB. 

With the unitary chiral rotation ip — > exp{— i'j 5 ^/ (2F)}iIj [24j-j2H| the Lagrangi- 
an (J2J) transforms into a Weinberg-type form C w containing the axial-vector coupling 



U+W eU-W 

« + 7 5 = 
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and the Weinberg- Tomozawa term [TT] : 

C w (x) = C {x) + Cf(x) + o(tt 2 ), (4) 
Co(x) = 4>(x){i - S(r) - 7 V(r)}^(x) - ^iz(x) (□ + M 2 ) n(x), 

£f{x) = -^pd^(x)i){x)^^Ti)(x) - |||7r i (x)^7r i (x)^(x)7 M r fc ^(a;), 

where Cj(x) is the 0(tz 2 ) strong interaction Lagrangian and □ = <9 M <9 M . 

In our calculation we do perturbation theory in the expansion parameter 1/F (where 
F is the pion leptonic decay constant in the chiral limit). We also treat the mass term of 
the current quarks as a perturbation. Dressing the baryonic three-quark core by a cloud 
of Goldstone mesons corresponds effectively to the inclusion of sea-quark contributions. 
All calculations are performed at one loop or at order of accuracy 0(1/ F 2 , in). 

We expand the quark field ip m the basis of potential eigenstates as 

i/;(x) = Yl b a u a (x)e- i£at + d\v p {x)^\ (5) 

a (3 

where the expansion coefficients b a and d\ are the corresponding single quark anni- 
hilation and antiquark creation operators. The set of quark {u a } and antiquark {vp} 
wave functions in orbits a and (3 are solutions of the static Dirac equation: 

[-i 7 ° 7 • V + 7 °S(r) + V(r) - £ a ]u a (x) = 0, (6) 

where £ a is the single quark energy. 

The unperturbed nucleon state is conventionally set up by the product of spin-flavor 
and color quark wave functions, where the nonrelativistic single quark wave function is 
replaced by the relativistic solution u (x) in the ground state. 

For a given form of the effective potential V e g(r) the Dirac equation in equation © 
can be solved numerically. Here, for the sake of simplicity, we use a variational Gaussian 
ansatz for the quark wave function given by the analytical form: 

u (x) = N exp ^ - ip l g_x \ XsXfXc (7) 

where N = [n 3 / 2 R 3 (l + f p 2 )] is a constant fixed by the normalization condition 
J d 3 x Uq(x) u (x) = 1 ; Xsi Xfi Xc are the spin, flavor and color quark wave functions, 
respectively. Our Gaussian ansatz contains two model parameters: the dimensional 
parameter R and the dimensionless parameter p. The parameter p can be related 
to the axial coupling constant calculated in zeroth-order (or three-quark core) 
approximation: 



2 



where 7 = (1 — |p )/(l + |p ) is a relativistic reduction factor. In our calculations 
we use a value of = 1.25 as obtained in the chiral limit of chiral perturbation 
theory The parameter R can be physically understood as the mean radius of the 
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three-quark core and is related to the charge radius of the proton in the leading-order 
(or zeroth-order) approximation as 



(rl) P LO = / d* X ul( X )x\ ( X ) = ^- tZ^L . ( 9 ) 



3R 2 1 + f p 2 
2 1 + |p 2 

The parameter R was fixed as R = 0.6 fm in our previous study of electromagnetic 
properties of nucleon which corresponds to (r%)^ = 0.6 fm 2 . 

The use of the Gaussian ansatz of equation (J7J) in its exact form restricts the scalar 
and the vector part of the potential to 

The above expressions are introduced as a static mean field potential confinement 
in PCQM, hence covariance cannot be fulfilled. As a consequence matrix elements are 
frame dependent: both Galilei invariance of the zeroth order baryon wave functions 
and Lorentz boost effects, when considering finite momenta transfers, are neglected. 
Approximate techniques j2Hl l2Hj have been developed to account for these deficiencies 
in static potential models. However, these techniques do not always agree and lead 
to further ambiguities in model evaluations. Furthermore, existing Galilean projection 
techniques are known to lead to conflict with chiral symmetry constrains JU]- I n the 
present manuscript we completely neglect the study of these additional model dependent 
effects. We focus on the role of meson loops, which, as shown in the content of cloudy 
bag model are not plagued by these additional uncertainties. 

According to the Gell-Mann and Low theorem we define the expectation value of 
an operator O in the PCQM by 

00 -n r r 

(6) = B (M J2 ~ / ^ ■■■ d A x n T[cY( Xl ) . . .cY W6]|^,>f • (12) 

n=0 U - J J 

The subscript c refers to contributions from connected graphs only and Cj(x) is the 
pion-quark interaction Lagrangian as already indicated in equation (j3J). The superscript 
B in equation (|T2*|) indicates that the matrix elements are projected on the respective 
baryon states. The projection of "one-body" diagrams on the nucleon state refers to 

xUU f ' fjS ' S XfXs ^ (N\J2(IJ)^\N), (13) 

i=i 

where the single-particle matrix element of the operators I and J, acting in flavor and 
spin space, is replaced by the one embedded in the nucleon state. For "two-body" 
diagrams with two independent quark indices i and j the projection prescription reads 
as 

xUl> A' 7 J!'" XfXs O xlxl I k 2 ' k 4 a XkX, ^ (N\ J> JO® ® (h J2) U) \N). (14) 



Axial form factor of the nucleon in the perturbative chiral quark model 



7 



We evaluate equation (fT2j) using Wick's theorem and the appropriate propagators. 
For the quark field we use a Feynman propagator for a fermion in a binding potential. 
The quark propagator iG^(x,y) is given by 

iG^(x,y) = ((j>o\T{ip(x)$(y)}\(j) ) -> ^2u a (x)u a (y) exp[-i£ a (x -y )]9(x -y Q ), (15) 

a 

where we restrict to the quark states propagating forward in time. The explicit form of 
the excited quark wave functions are obtained analytically as given in Appendix. For 
the mesons we use the free Feynman propagator for a boson field with 

iA^x -y) = (OlT^x^llO) = 5 l3 J ^Z^-^ZfJ - ( 16 ) 

To redefine our perturbation series up to a given order in terms of renormalized 
quantities a set of counterterms 5C has to be introduced in the Lagrangian. Thereby, 
the counterterms play a dual role: (i) to maintain the proper definition of physical 
parameters, such as nucleon mass and, in particular, the nucleon charge, and (ii) to 
effectively reduce the number of Feynman diagrams to be evaluated. Here we follow the 
formalism set out in reference [2J, but where in the present work intermediate excited 
quark states are included in the loop diagrams. In the following we attach the index 
"0 " to the renormalization constants when we restrict to the contribution of the ground 
state quark propagator and the superscript " F " when excited quark states are included. 

First we introduce the renormalized quark field ip r (x). It can be expanded in a set 
of potential eigenstates which are solutions of the renormalized Dirac equation with the 
full renormalized quark mass of 

3/1 \ 2 v- f°° „ , 2 1 



in). = in ( ^ / dk k 

1 ' UttFJ 



7 V47tF/ ^J co(k 2 )(uj(k 2 ) + AS a ) 

^(^[(fcl^^VjJ^iW+^l^f.J^Ji)] . (17) 

The expression for the renormalized quark mass includes self-energy corrections of the 
pion cloud, where 



F Ia (k) =N N a k 



drr 2 (g (r)g a (r)-f (r)f a (r)) / dfte 1 ™ C a Y la0 

Jn 

8 f°° f 
-2i— drr(/ (r)/«(r)) ^ dQ cosd e ikrcosd C a Y la0 (9, 4 

F IIa (k) = N N a -^ J™ drr(g (r)f a (r) - f (r)g a (r)) J d{le ikrcosd C a Y la0 (9, 0), (19) 

with the pion energy u(k 2 ) = a/ M 2 + k 2 ; k — \k\ is the pion momentum and A£ a = 
£a~£o is the excess of the energy of the quark in state a with respect to the ground state. 
The label a = (nl a jm) characterizes the quark state (principal quantum number n, non- 
relativistic orbital angular momentum l a , total angular momentum and projection j,m). 
For the Clebsch-Gordan coefficients we use the notation C a = {l a ^\\\j\) and Yi aQ (9, <f>) 
is the usual spherical harmonic. The explicit form of the radial wave functions g a (r) 



Axial form factor of the nucleon in the perturbative chiral quark model 8 

and f a {r), of the normalization constants (N a ) and of the energy difference (AS a ) are 
given in Appendix. 

When restricting the quark propagator to the ground state the expression above 
for the renormalized quark mass reduces to 

*5=*-^©'r d ** 4 w< <2o) 

where F n NN{k 2 ) is the 7rNN form factor normalized to unity at zero recoil (k 2 = 0): 



F 1TNN (k 2 ) = exp 



k 2 R 2 



k 2 R 2 



1 



(21) 



4 

In a second step we renormalize the effective Lagrangian including a set of 
counterterms. The renormalized interaction Lagrangian C^ r = C^.^ str + £^ r ' em contains 
a part due to the strong interaction, 

c W;str = c W;str + 5£ W;str^ (22) 

and a piece due to the electromagnetic interaction, 
The strong interaction term £}^> str i s given by 

The interaction of pions and quarks with the electromagnetic field is described by ^1] 



vr,(x)(9%-(x) - ^M^ {x)l ^ T . r{x) 



(25) 



which is generated by minimal substitution with 

d^ r - D^ r = d^ r + \eQA™ij> r , - D piX{ = + e^A^-, (26) 

where Q is the quark charge matrix. The set of counterterms, denoted by SC w ' ,str and 
fi£W;em are explained and given in reference j^]. 

Now we consider the nucleon charge and prove that the properly introduced 
counterterms guarantee charge conservation. Using Noether's theorem we first derive 
from the renormalized Lagrangian the electromagnetic current operator: 

+ + (27) 
It contains the quark component (j^ r ), the charged pion component the quark-pion 
component (j^ r7r ) and the contribution of the counterterm (6j!t r ): 

fir = VYQr = I (2flW - dr^dr) , 

# = E aij ir i & i ir j = 7r"i^7r + - Tr+i^Tr", 

= (^3 " 7TT7T ) ^ - ew^-f-fTdf', 

5j;=f(Z-l)/Qf 
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Here, Z is the renormalization constant determined by the nucleon charge conservation 
condition In the one- loop approximation following diagrams contribute to the 

nucleon charge [Figs.la-lf]: the three-quark diagram [Fig.la] with an insertion of the 
quark current j^ r , the three-quark diagram [Fig. lb] with the counterterm 5j!^ r (three- 
quark counterterm diagram), the self-energy [Figs.lc and Id], the vertex correction 
diagram [Fig.le] with the quark current j!t r , and finally the meson-cloud diagram [Fig. If] 
generated by the pion current j%. 

The analytical expression for the full renormalization constant, Z F is 

yF 1 _ 3 V [°° dk k 2 - 

( 4vFj p)2 J o u u}{k 2 ){u{k 2 ) + AS a ) 2 

x UjA;)F;jA;)-2^(fc 2 )F / Jfc)F; / JA ; )+^(A; 2 )F /7 Jfc)F; 7Q (fc)l . (29) 



When restricting intermediate quark states to the ground state equation (|2*9*|) yields the 
result 

= C?SY (30) 



400 \nF J J u 3 (k 2 ) 

We obtain a value of Z° = 0.9 jH] for our set of parameters. Inclusion of the excited 
quark states changes the value of the renormalization constant to a value of Z F = 0.7. 

3. The axial form factor of nucleon 

For the present purposes we have to construct the partially conserved axial-vector 
current A?: 

At = Fd^ + rn 5 ^ r - ^ r Yr^ r n k + (tt t tt, - tt 2 n ) ^ (31) 

+fr(Z - l)7 M 7 5 ^ r + o(n 2 ) . 

The axial form factor Ga{Q 2 ) of the nucleon is defined by the matrix element of the 
i — 3 isospin component and the spatial part of the axial vector current evaluated for 
nucleon states. In the Breit frame Ga{Q 2 ) is set up as [TH] : 

N s , (|) | J dW^aOjiV; (-|)\ = x \ a a N T l XNs G A {Q 2 ), (32) 

with Q 2 = —q 2 . Here, xn s and xn s , are the nucleon spin wave functions in the initial 
and final states; cr N is the nucleon spin matrix and rjy is the third component of the 
isospin matrix of the nucleon. At zero recoil (Q 2 = 0) the axial form factor satisfies the 
condition GU(0) = £U, where qa is axial charge of the nucleon. 

In the PCQM the axial form factor of the nucleon up to one loop corrections is 
initially given by 

.3 



n=0 ' ^ 

xT[CY' str (x 1 ) . . . £f' str (x n ) A 3 (x)} |0 o )f , (33) 
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with the interaction term 



cY ]Str (x 



d^{x)^{x)n 5 r^{x) - ^n^d^x^ix^n^ix) 



(34) 



where the superscript r refers to renormalized quantities. For the final calculation we 
include a set of excited states up to 2fuj in the quark propagator: the first p-states (lpi/2 
and lpz/2 in the non-relativistic spectroscopic notation) and the second excited states 
(1^3/2, 1g?5/2 and 2si/ 2 ). In other words we include the excited states whose energies 
satisfy to the restriction £ < A = 1 GeV (see discussion in the Sect.l). The diagrams 
to be evaluated are shown in Fig.2. 

Next, we present the analytical expressions for the axial form factor of the nucleon 
obtained in the PCQM. We start with the simplest case, where the quark propagator 
is restricted to the ground state contribution. The axial form factor of the nucleon is 
a sum of terms arising from different diagrams: the three-quarks diagram [Fig. 2a], the 
counterterm [Fig. 2b], the self-energy diagrams [Figs. 2c and 2d], the exchange diagram 
[Fig.2e] and the vertex-correction diagram [Fig.2f]. Other possible diagrams at one loop 
are compensated by the counterterm. The corresponding analytical expressions for the 
relevant diagrams are given in the following. 

(a) For the three-quark diagram (3g) [Fig. 2a] we obtain: 



G A (Q' 



3q 



LO 



Ga(Q ) + G A (Q 



3q 



NLO 



3q 



(35) 



where Ga^Q 2 )^ is the leading-order term (LO) evaluated with the unperturbed quark 



wave function (w.f.) 



u (x); G A{Q 2 )\^ q LO is a correction due to the renormalization of 
the quark w.f. uq(x) — > u r Q (x\rrf) which is referred to as next-to-leading order (NLO): 

LO 



g a (q 2 : 
g a (q 2 : 



9a Fttnn(Q 2 



3q 

NLO 
3q 



(36) 



pR 



1 9 2 

1+ 2 P 



G A {Q* 



LO 



3q 



(1 + fP 2 )' 

12 (2 - 3p 2 ) -4(1 + hp 2 ) Q 2 R 2 + p 2 Q 4 R 4 

The modified quark w.f. Uq(x; rh r ) is given by jH] 

Uq{x) + 5u r (x; m r ), 



5 

72 



exp 



Q 2 R 2 



where 



u {x; m 



6uq{x] rh r 



m 



pR 



2 l + |p 2 



1 + 21.2 

2 4' 



x c 

V 7 

R 2 1 



u (x 



(b) The three-quark counterterm (CT) [Fig. 2b] results in the expression: 

LO 



G A (Q' 



CT 



3q 



(37) 



(3f 



(39) 



(40) 



(c) The self-energy diagram I (SE;I) [Fig. 2c] yields 

SE-I 



Ga(Q 2 



pR 



(2 + 3p 2 ) \2tvF 



1 



dA; r 



(k 2 



o 



■D(fc,Q 2 ),(41) 
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V(k,Q 2 ) = exp 



3/2 



X 



(k + ^) 2 R 2 

k 2 Q 2 R 4 

2 + k^R 2 + exp (k^CpR 2 } (-2 + k^R 2 ^ 



(d) For the self-energy diagram II (SE;II) [Fig. 2d] we also obtain: 



G A {Q 2 ) 



SE-II 



G A {Q 2 ) 



SE-I 



(e) For the exchange diagram (EX) [Fig.2e] we get: 



G A {Q 2 ) 



_96 (o) pR f 1 
ex 5 9a (2 + 3p 2 ) \2irF 



dkk 



4 " nNN 



(k 2 ) 



L0 2 (k 2 ) 

(f) The vertex-correction diagram (VC) [Fig.2f] gives the contribution: 



(42) 
(43) 

V(k,Q 2 ).(U) 



G A {Q 2 ) 



vc 



3 > g m 



100 



1 



2ttF 



ttNN 



(Q 2 ) 



F 2 (k 2 ) 



cu 3 (F) 



• (45) 



In the next step we extend the formalism by also including excited states in the 
quark propagator. The leading-order expression of the three-quark diagram and the 
exchange term remain the same. In turn the following contributions must be extended. 

(a) In the three-quark NLO expression the appropriate renormalized mass has to 
be inserted with 



G A (Q 2 



NLO 



3q 



3 A 



■to, 



pR 



2 r (l + lp 2 ) 
5 



q \ lo 

i+y )g a (q 2 ) 

Z / 3<j 



72 



12 (2 - 3p 2 ) -4(1 + 5p 2 ) Q 2 R 2 + p 2 Q 4 R 4 



exp 



Q 2 R 2 



(46) 



(b) For the three-quark counterterm (CT) the renormalization constant has to be 
replaced accordingly 

LO 



G A (Q ) 



CT 



(Z* - 1)G A (Q 2 



3q 



(47) 



(c) For the self-energy diagram I (SE;I) we obtain the full expression 

2 



G A (Q 2 



SE;I 



10 

~3~ 



4ttF 



~d«^ ( * 2)f, "- w " f,, " ( * ) 



x 



dx 



where 

F IIIa (k_) = N N o 



_d_ 



oj(k 2 )(u(k 2 ) + AS a ) 
2k(l - x 2 )F IIIa (k_) + (yWx + (1 - 2x 2 )A;) F 7Va (£;_) 



(48) 



drr (g {r)f a (r)) dQe 



Ak— r cos ( 



F jya (£;_) = N N a J drr(f (r)g a (r)-g (r)f a (r)) 



J 

Jq 



x / dfie 



- ./o 

ifc_ r cos 6 



C a Y laO (0A), 



k 2 ± = k 2 + Q 2 ± 2k^/Q 2 x. 



(49) 
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(d) For the self-energy diagram II (SE;II) we get 

2 



G A (Q 2 



SE-JI 



10 

y 



i 



4ttF 



E 



™ ikk Mk 2 )F},Sk)-Fljk) 







X 



1 dx 2k(l-x*)F Va (k, 

-i 



oo(k 2 )(uj(k 2 ) + AS a ) 
- (^x + k)F IVa (k, 



(50) 



where 

F Va (k + ) = N N a 



_d_ 

dk 4 



drr(f (r)g a (r)) / d{le ik + r cose C a Y la0 (9, 0). 



(51) 



(e) For the vertex-correction diagram (VC) inclusion of excited states results in 



G A (Q 2 ) 



5^(Q 2 ) 



vc 



T- 

X 



dA; k 2 



+ 



uj(k 2 ){u{k 2 ) + A£ a ) (u(k 2 ) + A£ (3 )_ 
[u\k 2 ) (F IIa {k)F} h {k))-u{k 2 ) (F Ha (k)Fl(k) + F Ia (k)F} l0 (k)) 
(Fr a (k)Fl(k) 



where 



FaAQ 2 ) = ^^l™ drr 2 (^ i/3 (r) +2^ i/3 (r)), 

Aa,p(r) = {g a {r)gp{r) - f a {r)fp{r)) I dfi exp (i a/QV cos 9 j C a/3;1 ( 

B a , p {r) = f a (r)f p {r) J dQ exp (iy/Q*r costf) 



(52) 

(53) 
(54) 



x 



cos 2 9 C a p ; i{6, 0) + sin 9 cos 9 C a(3 . 2 (9, 0) 



(55) 
(56) 
(57) 



C Q/3;1 (0, 0) = C Q C^ aO (0, 0)^ o (^, 0) - DvDpY^O, <i>)Y lpl {0, 0), 

C a/3;2 (#, 0) = C Q D^ aO (0, 0)V(0> ^) e_i0 + DaCpY^O, <P)Y h0 (9, 0)e^, 

where _D a = (Z a l| — L an d //3 are the orbital quantum numbers of the intermediate 

states a and f3, respectively. 

The Q 2 -dependence (up to 0.4 GeV 2 ) of the axial form factor of the nucleon 
are shown in Figs. 3, 4 and 5, the description of each figure is given below. Due 
to the lack of covariance, the form factor can be expected to be reasonable up to 
Q 2 < p 2 = 0.4 GeV 2 , where p is the typical three- momentum transfer which defines the 
region where relativistic effect < 10% or where the following inequality p 2 / (4m%) < 0.1 
is fulfilled. 

The first result for the Q 2 -dependence of the axial form factor G A (Q 2 ) of the nucleon 
is indicated in Fig. 3. The numerical values are obtained, when truncating the quark 
propagator to the ground state or equivalently to the intermediate nucleon and delta 
baryon states in loop diagrams. Thereby, we also give the individual contributions of 
the different diagrams of Fig.2, which add up coherently. The leading order three-quark 
diagram dominates the result for the axial form factor, whereas pion cloud corrections 
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add about 20% of the total result. Here, both the exchange and self-energy terms give 
the largest, positive contribution. 

In a next step we include the intermediate excited quark states with quantum 
numbers lpi/2, J-P3/2, 1^3/2, 1^5/2 an d 2sy 2 in the propagator. The resulting effect on 
Ga{Q 2 ) is given in Fig. 4. We explicitly indicate the additional terms, which are solely 
due to the contribution of these excited states. The previous result, where the quark 
propagator is restricted to the ground state, is contained in the curve denoted by 
Total(GS). The inclusion of the intermediate excited states tends to induce a cancellation 
of the original pion cloud corrections generated for the case of the ground state quark 
propagator, thereby regaining approximately the tree level result. In Fig. 5 we give for 
completeness the full result for Ga(Q 2 ) including excited states in comparison with 
experimental data and with the dipole fit using an axial mass of Ma = 1.069 GeV and 
normalized to Ga{0) = 1-267 at zero recoil. The model clearly underestimates the finite 
Q 2 -behavior, but it should be noted that a similar effect occurs in the discussion of the 
electromagnetic form factors of the nucleon JO,; . The stiffness of the form factors can 
be traced to the Gaussian ansatz of the single quark wave functions and can also be 
improved when in addition resorting to a fully covariant description of the valence quark 
content of the nucleon [3U]. Hence the applicability of the PCQM is mostly for static 
quantities and low Q 2 observables of baryons. 

For the comparison with the data near Q 2 = we first turn to the results for the 
axial charge, Qa- In Table 1 we list the numerical values for the complete set of Feynman 
diagrams (Fig. 2), again indicating separately the contributions of ground and excited 
states in the quark propagator. 

The prediction for the axial charge including loop corrections are relevant for several 
reasons: 1) the tree level result for qa was previously adjusted to fix one (p) of the 
parameters. Since loop corrections essentially do not change this results, the previous 
model predictions remain meaningful; 2) the predicted small loop corrections to gA are 
consistent with similar results in chiral perturbation theory; 3) the generic role of excited 
states in loop diagrams are rather relevant in understanding the nucleon properties. This 
role was already exemplified in the case of N — A transition and sigma-terms 
and again is demonstrated in the case of gA- 

For low- momentum transfers, that is Q 2 < 1 GeV 2 , the axial form factor can be 
represented by a dipole fit 

Ga(Q 2 ) = ~^Li (58) 

in terms of one adjustable parameter Ma, the axial mass (or sometimes dipole mass). 
Therefore, the axial radius can be expressed in terms of the axial mass with: 



1 dG A {Q 



2 \ 



12 



{Ta) ~ 6 G A (0) dQ 2 Q2 =o ~ M 2 A m 
A comparison of the experimentally deduced values for the axial mass and the axial 
radius with our model results is given in Table 2. 
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4. Summary and Conclusions 

In summary, we have evaluated the axial form factor of the nucleon and, more important, 
its low Q 2 limits, such as the axial charge and the axial radius using a perturbative chiral 
quark model as based on an effective chiral Lagrangian. Since the PCQM is a static 
model, Lorentz covariance cannot be fulfilled. Approximate techniques to account for 
Galilei invariance and Lorentz boost effects were shown to change the tree level results 
by about 10%. Higher order, that is loop contributions, are less sensitive to these 
correction. The derived quantities contain, in consistency with previous works, only 
one model parameter R, which is related to the radius of the three-quark core, and are 
otherwise expressed in terms of fundamental parameters of low energy hadron physics: 
weak pion decay constant, and set of QCD parameters. In addition, another parameter 
(p of equation (6)), which is related to the amplitude of the small component of the 
single quark wave function, was originally set up to reproduce the value for the axial 
charge in the chiral limit with = 1-25 27\. Predictions are given for the fixed values 
of model parameters p and R in consistency with previous investigations. In particular, 
our result for the axial charge, gA = 1.19, is in reasonable agreement with the central 
value of data: g& = 1.267 ± 0.003. Thereby, contributions of excited quark states in 
loop diagrams play a considerable role in order to generate a small correction to the tree 
level result, which is required to account for the data point. This result, obtained in the 
context of the PCQM, is rather encouraging. Minor pion cloud corrections to the tree 
level result of gA justify in turn the appropriate choice for p or for used in previous 
works. Also, recent calculations of the axial charge up to order p A in chiral perturbation 
theory [3*1*1 l3*2*j imply rather large p 3 corrections leading to rather large uncertainties 
when going to the next order in the chiral expansion. Our model result can naturally 
explain the small correction to the one obtained in the chiral limit, but only when going 
beyond nucleon and delta states in the loop diagrams. 

Acknowledgments 

This work was supported by the Deutsche Forschungsgemeinschaft (DFG Nos.FA67/25- 
3, GRK 683). K. K., S. C. and Y. Y. acknowledge the support of Thailand Research 
Fund (TRF, Grant No.RGJ PHD/00187/2541, PHD/00165/2541) and the Deutscher 
Akademischer Austauschdienst (DAAD, Grant No. A/01/19908). K. P. thanks 
the Development and Promotion of Science and Technology Talent Project (DPST), 
Thailand for financial support. 

Appendix: Solution of Dirac equation for the effective potential 

In this appendix we indicate the solutions to the Dirac equation with the effective 
potential V e g(r) = S(r) + 7 c V(r) , with r = \x\. The particular choice of scalar S(r) 
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and the time-like vector V(r) parts are given by 

l-3p 2 



S(r) 



V(r) = S 



2pR 2R 3 
l + 3p 2 



P 2 

r . 
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(60) 
(61) 



2pR 2R? 

The quark wave function u a (x) in state a and eigenenergy £ a with the specific choice 
of Ves satisfies the Dirac equation 

[-i 7 °7 • V + 7 °5(r) + V(r) - £ a ]u a {x) = 0. (62) 

The solution of the Dirac spinor u a (x) to (Jfi2*|) can be written in the analytical form |17j : 



u n [x) 



9a(r) 
icr ■ ff a (r) 



y a (r)Xf Xc 



(63) 



For the particular choice of potential the radial functions g a (r) and f a (r) satisfy the 
form 



9a{r) = 
where for the j = I + \ 

fair) = pa 

and for j — I — ~ 



r 

Rn 



I / 2 

T i + i/2 ( r l 

Ln - 1 [ K 



fair) 



-Pc 



r 

R a 



i-i 



r 



n + / 



i+i 



■ 1+3/2 



R 2 7 



L 



1-1/2 
n-1 



k) 



\R 2 



_ r 2 



(64) 



(65) 



(66) 



The label a = (nljrrij) characterizes the state with principle quantum number n = 
1, 2, 3, orbital angular momentum I, total angular momentum j = Z±| and projection 
rrij. Due to the quadratic nature of the potential the radial wave functions contain the 
associated Laguerre polynomials L^(x) with 

in + k)\ 



m=0 



(n — m)\{k + m)!m! 
The angular dependence, y a ir) = y^m^ir), is defined by 

yijm s {r)= ^2 \ lrni o m *i m 3 / Y imi(r)xh 



mi,m s 



(67) 



(6J 



where Y; m; (f) is the usual spherical harmonic. Flavor and color part of the Dirac spinor 
are represented by Xf an d Xc, respectively. 

The two coefficients R a and p a in the redial function of state a are of the form 

R a = R(l + AS a pRy lf \ (69) 



Pa = P 



Ra 



(70) 
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and are related to the Gaussian parameter p, R of equation (J7J) . The quantity A£ a = 
£ a — £q is the excess of the energy of the quark state a with respect to the ground state. 
A£ a depends on the quantum number n and / and is related to the parameters p and 
Rhy 

2 



A£ a + 



3p 
R 



AS a + 



1 
~pR 



_P_ 

R 3 



(in +21 - If 



The normalization constant, which results in 



2 -2(n + ^i/2)^i/2^3 (2n + 2Z)! L + #( 2 



N 

'(// • /)!(;/ /)! 
is obtained from the normalization condition 



Pi[2n + l- 



-1/2 



t>oo 

I d^XU^ix^ai^x) = 1. 

Jo 



(71) 



(72) 



(73) 
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Table 1. Contributions of the individual diagrams of Fig. 2 to the axial charge qa- 
Separate results for the inclusion of ground (GS) and excited states (ES) in the quark 



propagator are indicated. 


9A 


GS quark propagator 




3q-core 




LO 


1.25 


NLO 


-0.06 


Counterterm 


-0.12 


Exchange 


0.23 


Vertex correction 


0.01 


Self-energy 


0.19 


GS contribution 


1.50 


ES quark propagator 




NLO 


-0.31 


Counterterm 


-0.25 


Vertex correction 


0.03 


Self-energy 


0.22 


ES contribution 


-0.31 


Total(GS+ES) 


1.19 


Experiment 2 


1.267 ±0.003 



Table 2. Comparison of the axial mass and the axial radius between experimental 
values and the result from the PCQM. 





Model 


Experiment 


M A (GeV) 


0.78 


1.069 ±0.016[T] 


(r A ) 1/2 (fm) 


0.88 


0.639 ±0.010[T] 
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Figure 1. Diagrams contributing to the nucleon charge: triangle diagram (a), 
triangle counterterm diagram (b), self-energy diagrams (c) and (d), vertex correction 
diagram (e) and meson-cloud diagram (f). 
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Figure 2. Diagrams contributing to the axial form factor of the nucleon : three-quark 
core (a), counterterm (b), self-energy (c and d), exchange (e) and vertex correction (f). 
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Figure 3. Model result for the axial form factor of the nucleon Ga(Q 2 )- The coherent 
contributions of the different diagrams of Fig. 2 are indicated when restricting to the 
ground state(GS) quark propagator. 
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3q-NLO(ES) 

CT(ES) 

VC(ES) 

SE(ES) 

— — Total(ES) 
Total(GS+ES) 



-0.25 



-0.50 




Figure 4. Model result for the axial form factor of the nucleon Ga(Q 2 ) when excited 
states are included in the quark propagator. The full ground state result is contained 
in the curve labelled by Total(GS). Excited state (ES) contributions of the individual 
diagrams are indicated separately. 




Figure 5. The axial form factor of the nucleon in the PCQM in comparison with a 
dipole fit (axial mass Ma = 1.069 GeV) and with experimental data. Data are taken 
from references |55]-|55). 



